Search for: All records

In Paper I [H. Francisco, A. C. Cancio, and S. B. Trickey, J. Chem. Phys. 159, 214102 (2023)], we gave a regularization of the Tao–Mo exchange functional that removes the order-of-limits problem in the original Tao–Mo form and also eliminates the unphysical behavior introduced by an earlier regularization while essentially preserving compliance with the second-order gradient expansion. The resulting simplified, regularized (sregTM) functional delivers performance on standard molecular and solid state test sets equal to that of the earlier revised, regularized Tao–Mo functional. Here, we address de-orbitalization of that new sregTM into a pure density functional. We summarize the failures of the Mejía-Rodríguez and Trickey de-orbitalization strategy [Phys. Rev. A 96, 052512 (2017)] when used with both versions. We discuss how those failures apparently arise in the so-called z′ indicator function and in substitutes for the reduced density Laplacian in the parent functionals. Then, we show that the sregTM functional can be de-orbitalized somewhat well with a rather peculiarly parameterized version of the previously used deorbitalizer. We discuss, briefly, a de-orbitalization that works in the sense of reproducing error patterns but that apparently succeeds by cancelation of major qualitative errors associated with the de-orbitalized indicator functions α and z, hence, is not recommended. We suggest that the same issue underlies the earlier finding of comparatively mediocre performance of the de-orbitalized Tao–Perdew–Staroverov–Scuseri functional. Our work demonstrates that the intricacy of such two-indicator functionals magnifies the errors introduced by the Mejía-Rodríguez and Trickey de-orbitalization approach in ways that are extremely difficult to analyze and correct.

 

The revised, regularized Tao–Mo (rregTM) exchange-correlation density functional approximation (DFA) [A. Patra, S. Jana, and P. Samal, J. Chem. Phys. 153, 184112 (2020) and Jana et al., J. Chem. Phys. 155, 024103 (2021)] resolves the order-of-limits problem in the original TM formulation while preserving its valuable essential behaviors. Those include performance on standard thermochemistry and solid data sets that is competitive with that of the most widely explored meta-generalized-gradient-approximation DFAs (SCAN and r2SCAN) while also providing superior performance on elemental solid magnetization. Puzzlingly however, rregTM proved to be intractable for de-orbitalization via the approach of Mejía-Rodríguez and Trickey [Phys. Rev. A 96, 052512 (2017)]. We report investigation that leads to diagnosis of how the regularization in rregTM of the z indicator functions (z = the ratio of the von-Weizsäcker and Kohn–Sham kinetic energy densities) leads to non-physical behavior. We propose a simpler regularization that eliminates those oddities and that can be calibrated to reproduce the good error patterns of rregTM. We denote this version as simplified, regularized Tao–Mo, sregTM. We also show that it is unnecessary to use rregTM correlation with sregTM exchange: Perdew–Burke–Ernzerhof correlation is sufficient. The subsequent paper shows how sregTM enables some progress on de-orbitalization.

 

Abstract The hybrid design of the Pierre Auger Observatory allows for the measurement of the properties of extensive air showers initiated by ultra-high energy cosmic rays with unprecedented precision. By using an array of prototype underground muon detectors, we have performed the first direct measurement, by the Auger Collaboration, of the muon content of air showers between $$2\times 10^{17}$$ 2 × 10 17 and $$2\times 10^{18}$$ 2 × 10 18 eV. We have studied the energy evolution of the attenuation-corrected muon density, and compared it to predictions from air shower simulations. The observed densities are found to be larger than those predicted by models. We quantify this discrepancy by combining the measurements from the muon detector with those from the Auger fluorescence detector at $$10^{{17.5}}\, {\mathrm{eV}} $$ 10 17.5 eV and $$10^{{18}}\, {\mathrm{eV}} $$ 10 18 eV . We find that, for the models to explain the data, an increase in the muon density of $$38\%$$ 38 % $$\pm 4\% (12\%)$$ ± 4 % ( 12 % ) $$\pm {}^{21\%}_{18\%}$$ ± 18 % 21 % for EPOS-LHC , and of $$50\% (53\%)$$ 50 % ( 53 % ) $$\pm 4\% (13\%)$$ ± 4 % ( 13 % ) $$\pm {}^{23\%}_{20\%}$$ ± 20 % 23 % for QGSJetII-04 , is respectively needed.